Synchrotron dynamics in Compton x-ray ring with nonlinear momentum compaction 
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The longitudinal dynamics of electron bunches with a large energy spread circulating in the storage 
rings with a small momentum compaction factor is considered. Also the structure of the longitudinal 
phase space is considered as well as its modification due to changes in the ring parameters. The 
response of an equilibrium area upon changes of the nonlinear momentum compaction factor is 
presented. 
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I. INTRODUCTION 



Engagement of electron storage rings for production 
of x rays through Compton scattering of laser photons 
against ultrarelativistic electrons was proposed in 1998 
PJ . Two basic schemes exist so far. One of them supposes 
use of electron beams with unsteady parameters [2( and 
applies the continual injection (and ejection of circulating 
bunches by the next injecting pulse) of dense intensive 
bunches. The second scheme is based on the continuous 
circulation of bunches. To make the bunches acquired a 
sufficiently large energy spread confined (see [3J), a lattice 
with a small controllable momentum compaction factor 
is proposed to employ Q. The longitudinal dynamics 
in the small compaction lattice is governed not only by 
the linear effects of the momentum deviation but by the 
nonlinear ones as well. 

In Compton sources storing the bunches with the large 
energy spread which can be as high as a few percents, 
ring's energy acceptance becomes compared to the en- 
ergy spread. To get proper life time of the circulating 
electrons, the energy acceptance a = ma,x(E — E s )/E s 
(E s is the energy of synchronous particle) should be high 
enough. 

Within a linear approximation according to the en- 
ergy deviation, the acceptance can be increased either 
by enhancement of the RF voltage, V r i, or by decreas- 
ing of the linear momentum compaction factor cto since 

O OC y/Vrf/oLQ. 

The paper presents results of study on the longitudi- 
nal dynamics of electron bunches circulating in storage 
rings with a small linear momentum compaction factor 
ato- Structure of the phase space are considered and its 
deformation with changes in the ring lattice parameters. 
In particular, the size of stable area as a function of the 
RF voltage and momentum compaction is evaluated. 



II. FINITE-DIFFERENCE MODEL OF 
LONGITUDINAL MOTION 

Let us consider a model of the ring comprised only two 
components: drift and radio frequency (rf) cavity. For 
the sake of simplicity we will suggest the cavity infinitely 
short, in which the particle momentum (energy) suffer 
an abrupt change while the phase of a particle remains 
unchanged. On the contrary, the phase of a particle trav- 
eling along the drift changes while the energy remains 
invariable. The longitudinal motion in a such idealized 
ring will be described in canonically conjugated variables 
4> (the phase about zero voltage in the cavity) and the 
momentum p = (7 — ^ s )lls equal to the relative devia- 
tion of the particle energy from the synchronous one (7 S 
is the Lorentz factor of the synchronous particle) . 

To study systems able to confine the beams with large 
energy spread, one needs to account not only the linear 
part of the orbit deviation from the synchronous one, but 
nonlinear terms as well: 



D x p + D 2 p 2 



(1) 



where D\ and D2 are the dispersion functions of the first 
and second orders, respectively. 

Accordingly, relative lengthening of a (flat) orbit is 



(2) 



where Lq is the length of synchronous orbit, p(s) the local 
radius of curvature, s the longitudinal coordinate. The 
coefficients cvq and &i arc determined as 
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In accordance with the definitions for ceo and ai, the 
momentum compaction factor a c can be written as 
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To study the phase dynamics in a storage ring with 
small momentum compaction factor a , the next terms 
of expansion of the compaction over the energy deviation 
should be accounted for, hence — the higher terms in the 
sliding factor »7 H, 0, • Magnitude of 77 characterizes 
a relative variation of the phase due to changes of the 
particle velocity and orbit length. It is determined by 
the relation 



— r = v(p)p ~ 0?o + viP + ■ ■ ■ )p j 



with ?7o and 771 having been determined by 

11 = a - I hi 

111 = ai + 770 
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(6a) 
(6b) 



The finite-difference equations for the phase <p and the 
variation of relative energy p in the model under consid- 
eration read 

<pf = <P% + {noPi + kiP^)At ; (7a) 
p f =p-i- U ri sin 0/ At , (7b) 



where 

(3c 

AT = Tf-n= —(tf -U) , 

the subscripts i and / correspond to the initial and 
final values, respectively. The dimensionless variable 
t = t(3c/L represents time expressed in number of rota- 
tions [t is time, (3c the velocity of a particle) . The factors 
kq and K\ at a large "f s are determined by the expressions 
Kq = 2nhrjQ w lixha^^ K\ = 27r/i77i ~ 2i:h(ao + an) ( h 
the harmonic number). 

From Eqs. Q, differential (smoothed) equations can 
be deduced. As it seen, the RHS of l|7b(l contains the 
final value of the phase cj> f expressed via the initial value 
<pi and momentum pi by the equation (|7a() . 

Let us expand sin^f into series of powers of At: 



sin^/ = sin [fa + (K Pi + kij>?)At] 

« sin^ + cos (ptinoPi + KiPi)Ar 



(8) 



Since At can not be regarded as infinitesimal (formally 
Eqs. present a complete turn, At = 1), then the 
linear term can be neglected if KoPi + K\pf <C 1. In the 
considered case it can be done since maximum of the 
energy spread does not exceed a few percents, and the 
momentum compaction factor ao supposed small. From 
these assumptions, finite difference equations reduce to 



Acp 
At" 



K0P1 + KlPi 



Ap . 

— = -£/ rf sm< 
At 



(9b) 



III. DIFFERENTIAL MODEL OF MOTION 

Noting of formal similarity of Eqs. JSJl to canonical 
Hamilton equations describing a mathematical pendu- 
lum, we can use a smoothed analog to these equations 
(a differential substitute for a finite difference equation, 
At — > 0) to facilitate analysis of the motion 



dT 

dp 
dr" 



= K p + KlP 

= —U r { sin <p 



(10a) 
(10b) 



A Hamilton function for (|10|) possesses a specific form 
with the cubic canonical momentum term 



H 



yP 3 + yP 2 + C/ lf (l-cos, 



(11) 



To analyze a phase portrait of the system, it is expe- 
dient to present Hamilton function of the longitudinal 
motion in the reduced form: 



H = f i 



P 3 



P 2 



1 — COS ( 



where 



k / 2TTha j s E 

TTP = \ 77 P ' 

t/rf V ev 

njUri _ (ao + ai) 2 eV 
Kq 2irhal^ s E 



(12) 

(13a) 
(13b) 



Phase portraits of motion with the Hamiltonian H12|) 
represented in Fig. ^ Topology of the phase plane is 
governed by the magnitude and sign of the parameter 
p. At zero value, [i — 0, the Hamiltonian (|llfl or 112|) 
has a form of mathematical pendulum; its phase plane is 
presented in Fig.^a). 

Within the interval < /i 2 < 1/12, there an additional 
area of finite motion appears; this area is separated from 
the main area with the band of infinite motion as depicted 
in Fig. Gib) . When the parameter /j, exceeds the critical 
value /j 2 = 1/12 [see Fig. [He)], e.g. 1/12 < /j 2 < 00, 
the structure of the phase plane will have changed as is 
represented in Fig. ^d). 

The dimension of a stable (finite) longitudinal motion, 
i.e., the area comprised by a separatrix, is in direct pro- 
portion with ratio of the ring parameters. For the con- 
sidered case of the nonlinear Hamiltonian (|12fl . the sep- 
aratrix height (size along the p axis) is determined by 
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FIG. 1: Phase portrait of longitudinal motion with account for the cubic nonlinearity at different values of the parameter u. 
(a): fj, = 0, (b): fi < |/x c |, (c): a = /i Cl (d): a > \fj, c \ 



for [i < n c <|14a[) and /i > fi c (|14b(l . respectively. 

The phase width of the separatrix (dimension along 
the 4> axis) is determined by expressions 



A(f> = 2vr ; 
A</> = 2 arccos 



7i7? 



3C^ r f (ao + ai) 



(15a) 
(15b) 



TABLE I: Ring parameters 



parameter 


desig 


value 


Accel, voltage (Volt) 


Kf 


4 x 10 5 


Lorentz factor 


Is 


84 


Harmonic number 


ll 


32 


Lin. comp. factor 




0.01 


Quad. comp. factor 


ai 


0.2 



for the subcritical and overcritical values of the parame- 
ter fJL. 

Dependence of the phase and momentum separatrix 
extensions on rf amplitude at fixed other parameters, 
which values are listed in Tab. [I] is presented in Fig. [21 

As it can be seen from the plot in Fig. [3 while in- 
creasing the parameter U r f the separatrix height grows 
up reaching its maximum, Ap k 7.1 X 10~ 2 , at f7 r f( c ) « 
3.8 x 10~ 4 (which is equal to the rf voltage of Vi c \ ~ 
16.3 kV at ls = 84). ' 



With further increase in the rf voltage, the separatrix 
height remains constant. 

The separatrix width remains constant with increase 
of the rf voltage up to the critical value t/ r f ( c ) , then it is 
diminishing. 

In Fig. |31 a dependence of the separatrix dimensions 
upon the linear momentum compaction factor under 
other system parameters fixed is presented. 

Quite the reverse to the dependence A0 = A0(/7 r f), 
a dependence of the separatrix width upon the linear 
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FIG. 2: Separatrix height (above) and width (below) as func- 
tions of the parameter U r t 



compaction factor, A<j) = A0(ao), is increasing while 
cto grows. At a certain critical value of the linear mo- 
mentum compaction factor a ( c ) (in the suggested case 
&o(c) ~ 0.03 ), the width of equilibrium area has reached 
its maximum and remains constant with further increase 
in «o. A dependence of the separatrix height on ao is 
of increasing within interval < «o < a ( c y Then, af- 
ter the maximum at ao — ao( c ) this dependence becomes 
declining, coming to zero at a large ao- 

Since the phase volume enclosed within the separatrix 
(and, therefore, the storage ring acceptance) is propor- 
tional to product of the transverse dimensions of the sep- 
aratrix, a ~ ApAcj), then from comparison of the plots in 
Fig- Eland Fig.|3]it follows that optimal working point is 
about the critical parameters. 

In addition, it can be seen that, dislike a linear lat- 
tice, nonlinear terms in the momentum compaction fac- 
tor restrict the infinite increase of energy acceptance 
with decreasing of the linear momentum compaction fac- 
tor: The acceptance increase takes place while the linear 
compaction is above certain critical value a ( c ), which is 
determined by the ring lattice parameters according to 



FIG. 3: Separatrix height (above) and width (below) as func- 
tions of Qo 

equality 

(ao + " l)2ey = -I . (16) 
2nha^ s Eo 12 v ; 

With further decrease of ao the acceptance decreases 
also. 

To validate use of differential (smoothed) equations of 
motion (fTU|) for analysis of Compton storage rings, a code 
has been developed based on the finite difference equa- 
tions 0. A simulated phase space portrait at the ring 
parameters listed in Tab. [I] for /i > /i c is presented in 
Fig.H 

From the figure it follows that the electrons can be 
confined within not only the "linear" area (minimum 
of Hamilton function (|12|) ). but the "nonlinear" as well. 
(The nonlinear stable region disappears in a linear lat- 
tice.) RMS sizes and the center of weight positions per- 
fectly correspond to the analytical estimations presented 
above. 



IV. SUMMARY. CONCLUSION 

Results of the study on dynamics of synchrotron mo- 
tion of particles in the storage rings with the nonlinear 
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FIG. 4: Distribution of confined electrons over the longitudi- 
nal phase plane in a system with cubic nonlinearity at /i > \i c ; 
left bunch corresponds to the "linear" case, right - to the 
"nonlinear" (additional) one 



momentum compaction factor presented in the paper, 
can be digested as follows: 

• Grounded on a simplified model of the storage ring, 
the finite-difference equations were derived. Hamil- 
tonian treatment of the phase space structure was 
performed. As was shown, the structure of the 
phase space is governed by ratios of the ring pa- 
rameters. An analytical expression for the factor 
which determines the topology of the longitudi- 
nal phase space, was derived. 



• Dependencies of the sizes of the equilibrium areas of 
the synchrotron motion in a nonlinear lattice were 
derived. Analysis of dependence of the longitudinal 
acceptance upon the amplitude of rf voltage, and 
the linear compaction factor at the fixed quadratic 
nonlinear term was presented. As was shown, the 
acceptance is growing up only to a definite magni- 
tude which determines by the critical value of pa- 
rameter (j, = fi c . It was emphasized that in order to 
maximize the acceptance of a lattice with a small 
linear momentum compaction factor and a wide en- 
ergy spread of electrons in the bunches, the system 
parameters should be chosen close to the critical 
value of \i. 

• To validate the use of smoothed equations of mo- 
tion, a simulating code was developed. The code 
is based on the finite-difference equations. The re- 
sults of simulation manifest a good agreement with 
the theoretical predictions on the sizes and position 
of equilibrium areas. 

The results obtained allow to make the following con- 
clusion: Enlargement of the energy acceptance of a ring 
by decreasing of the momentum compaction factor is lim- 
ited with the nonlinearity in the compaction factor. De- 
creasing of the linear compaction factor below the certain 
limit causes the reversed effect — decreasing of the ac- 
ceptance. 

Similar consequence corresponds to the build-up of the 
rf voltage: Increase of the voltage above a certain limit 
causes narrowing of possible bunch lengthes while the 
energy acceptance remains constant. This effect can lead 
to decrease in the injection efficiency for high rf voltages. 
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